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Abstract 


Few  studies  have  been  made  on  the  stress/strain  field  or 
plastic  zone  size  ahead  of  crack  tip  in  a  high  temperature 
environment  under  varying  load  frequencies  and  stress  levels. 
Fewer  studies  have  incorporated  compressive  loads  or  analyzed 
the  effect  of  a  negative  R-ratio  on  the  fatigue  characteristics 
of  the  superalloy  IN-100. 

This  study  involves  extending  existing  analysis  of  the 
stress  field  and  plastic  zone  ahead  of  a  crack  tip  in  a  com¬ 
pact  tension  specimen,  through  a  larger  number  of  load  cycles 
and  examining  the  nearly  unexplored  area  of  compressive  loading 
of  a  crack  in  a  superalloy.  A  USAF  Materials  Laboratory 
finite  element  computer  program  named  VISCO  was  used  for  this 

M.  - 

study.  The  Bodner-Partom  viscoplastic  constitutive  equations 
for  describing  the  material  behavior  were  utilized.  Load 
spectra  included  various  frequencies  with  R-ratios  of  0.1  and 


-1.0  (zero  mean  load) „ 


I  Introduction 


Current  emphasis  in  the  United  States  Air  Force  on  Product 
Assurance  and  equipment  reliability  coupled  with  the  ever 
present  need  for  cost  savings  has  created  the  problem  of 
maximizing  the  life  cycle  of  airframe  and  engine  components 
without  jeopardizing  mission  safety.  There  are  numerous 
standards  governing  the  design,  construction,  and  maintenance 
of  airframe  and  engine  components.  These  standards  include 
requirements  for  periodic  inspections  to  detect  flaws. 

In  contrast  to  the  normal  requirements  are  some  low  cycle 
fatigue-limited  jet  engine  parts  which  may  be  retired  from 
service  when  no  flaws  have  yet  been  found  in  them.  Jet  engine 
turbine  disks  are  removed  from  service  at  a  time  when  statis¬ 
tically  1  in  1000  would  be  expected  to  have  initiated  a  crack 
of  some  finite  length  (0.03  in).  Eighty  percent  of  the  disks 
have  at  least  10  useful  lifetimes  remaining,  though  no  attempt 
is  made  to  utilize  the  statistically  "failed"  disks  [1]„ 

From  a  mission  safety  standpoint,  this  retirement  approach 
has  been  very  successful.  This  approach  is  also  very  costly. 
It  has  been  estimated  that  replacement  costs  for  low  cycle 
fatigue-limited  jet  engine  components  are  in  the  vicinity  of 
$100,000,000  for  the  1980-1985  time  period  [2].  A  signifi¬ 
cant  reduction  of  this  cost  could  be  realized  if  a  procedure 
was  developed  to  utilize  a  disk  after  it  has  statistically 
"failed." 
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The  life  cycle  of  a  turbine  disk  is  a  complex  one  consist¬ 
ing  of  frequent  load  cycle  variations  with  ambient  temperatures 
of  up  to  1350  F.  The  elevated  temperatures  introduce  time 
dependent  creep  phenomena  which  interact  with  the  varied  load 
spectra  to  produce  complex  material  behavior.  If  the  material 
behavior  could  be  determined  for  the  typical  life  cycle  then 
accurate  remaining  life  predictions  could  be  made  for  compo¬ 
nents  with  subcritical  flaws.  Only  those  components  with  a 
quantifiable  critical  flaw  size  would  therefore  need  to  be 
retired  [3] . 

Approach 

This  study  attempts  to  quantify  the  material  behavior 
under  cyclic  loading  of  IN-100,  a  superplastically  forged 
nickel-based  superalloy  used  in  turbine  disks  for  the  F-100 
jet  engine.  VISCO,  a  computer  program  developed  by  Hinnerichs 
[4]  was  used  for  this  analysis.  VISCO  is  a  finite  element 
program  which  uses  constant  strain  triangular  elements  and 
has  the  capability  to  run  with  non-static  loads „ 

The  Bodner-Partom  flow  law  subroutine  in  VISCO  was  used 
to  model  the  plastic  flow  during  the  load  cycling.  The  Bodner 
material  parameters  for  IN-100  at  1350  F  were  experimentally 
determined  by  Stouffer  [5],  The  Bodner-Partom  viscoplastic 
flow  law  has  the  capability  to  predict  the  behavior  produced 
by  cyclic  effects,  time  dependent  creep,  strain  hardening. 


and  plastic  deformation.  This  flow  law  is  integrated  through 
time  by  an  Euler  extrapolation  scheme  and  incorporated  into 


the  finite  element  program  by  utilizing  the  residual  force 
technique  [4] . 

A  compact  tension  specimen  was  modeled  differently  than 
that  used  in  previous  work  [12] ,  to  allow  a  larger  number  of 
load  cycles  to  be  examined.  Load  was  input  as  a  sawtoothed 
stress-time  pattern  of  constant  maximum  and  minimum  stress. 
Load  ratios  of  0.1  and  -1.0  were  examined.  During  each  time 
step  the  stress,  strain,  plastic  work  and  the  Z  hardness 
material  parameter  were  calculated  for  each  element  in  the 
model. 

The  load  ratio,  herein  referred  to  as  R-ratio,  is  defined 
as  the  ratio  of  the  applied  minimum  load  to  maximum  load. 


Viscoplasticity  is  a  combination  of  two  strain  groups, 
time  dependent  nonrecoverable  strain  and  plasticity.  Time 


dependent  strains  accumulate  at  a  finite  rate  and  are  there¬ 
fore  time  dependent.  Plastic  strains  are  permanent  and 
develop  instantly  since  they  are  independent  of  time.  The 
combined  strain  effect  of  time  dependent  strain  and  plasticity 
is  measurable,  irreversible  and  can  be  utilized  in  a  unified 
plastic  strain  rate  model.  All  materials  exhibit  some  plastic 
behavior,  even  when  the  applied  loads  produce  stress  below 
the  yield  stress.  Most  plastic  deformation  below  the  yield 
stress  is  small  and  is  therefore  usually  neglected  [6] . 

Separating  the  total  strain  into  elastic  (reversible)  and 
plastic  (nonreversible)  components  yield: 


e  p 
e.  .  *  e .  .  +  e  . \ 
13  13  13 


(2.1) 


By  taking  the  time  derivative  of  equation  (2.1),  an  expres¬ 
sion  for  the  total  strain  rate  is  obtained. 


&.  .  *  e.®  +  e.? 

13  13  13 


(2.2) 


where  is  the  total  strain  rate,  is  the  elastic  strain 
rate  and  6^?  is  the  viscoplastic  strain  rate.  The  elastic 
strain  rate  is  related  to  the  stress  rate  by  the  time 

derivative  of  Hooke's  law.  The  viscoplastic  strain  rate 


assuming  incompressibility  and  isotrophy,  is  taken  to  follow 
the  Prandtl-Reuss  flow  law  of  classical  plasticity 


(2.3) 


where  are  the  components  of  the  deviatoric  stress  tensor 

and  X  is  a  scalar  that  reflects  the  viscosity  of  the  material. 
In  a  classical  theory,  plastic  deformation  begins  at  yield 
and  depends  on  a  yield  criterion.  It  can  be  shown  that,  at 
yield,  the  Prandtl-Reuss  relations  imply  the  Von  Mises  yield 
criterion  [7] „  The  Von  Mises  yield  criterion  is  used  primarily 
for  metals  behavior  and  is  based  on  the  distortion  strain 
energy  theory.  Yielding  begins  when  the  distortion  strain 
energy  equals  the  distortion  strain  energy  required  to  produce 
yielding  in  a  simple  uniaxial  stress  test. 

The  distortion  strain  energy  can  be  written  as 

U.  =  \S.  .de*.  .dVol  (2.4) 

d  1  i]  i] 

V 

where  de'^  and  are  tensors  related  only  to  distortion, 

(S. .  is  the  deviatoric  stress  tensor).  U.  can  be  written 
i]  d 

in  terms  of  the  second  invariant  of  deviatoric  stress,  J^r 
and  the  shear  modulus  G  as 
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I 

t 


where  J2  in  terms  of  the  principal  stresses  is 

J,  =  7  [  (a  -a  )2+(o  -0  )2+(o  -o  )2 

2  6  |_  1  2  2  3  IS 

The  distortion  strain  energy  can  now  be  written 

U.  =12G  f  (a  -a  )  2+  (a  -a  )2+(a  -a  )2 

a  |_  1  2  2  3  is 

Since  yielding  is  related  t\  the  strain  energy  in  a  uniaxial 
test,  and  can  be  reduced  accordingly.  For  a  unixial 

test 

o  =  a  a  =  a  =  0  (2.8) 

1  ys  23 

where  ayg  equals  the  stress  at  the  start  of  yielding.  J2 
reduces  to 


] 


(2.6) 
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(2.7) 
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a 

ys 


and  the  distortion  strain  energy  becomes 


(2.9) 


(2.10) 


Equating  the  multiaxial  distortion  strain  energy  to  that  of 
the  uniaxial  case  yields 


V 
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(2.11) 


where  k  is  a  proportionality  constant,  Multiaxial  plastic 
yielding  is  predicted  to  occur  when  J2  reaches  the  critical 
value  at  yield  in  a  uniaxial  stress  test. 


Bodner-Partom  Constitutive  Law 

Bodner  and  Partom's  constitutive  law  is  based  on  the  study 
of  dislocation  dynamics  [8],  The  constitutive  equations  are 
able  to  represent  the  principal  features  of  cyclic  loading 
behavior,  including  softening  upon  stress  reversal,  cyclic 
hardening  or  softening,  cyclic  saturation,  cyclic  relaxation, 
and  cyclic  creep.  The  formulation  of  the  equations  is 
arrived  at  by  squaring  the  Prandtl-Reuss  relation: 


(2.12) 


p 

where  D2  is  the  second  invariant  of  the  plastic  strain  rate 
and  J2  is  the  second  invariant  of  the  deviatoric  stress. 
The  plastic  deformation  rate  has  a  functional  relationship 
with  the  stress  invariants,  i.e. 


°2  =  D2  (J2'VT) 


(2.14) 
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Utilizing  the  above  equations,  the  parameter  X  can  be  found 
to  be: 


r  p  ~i  1/2 

X  -[d2  (J2,2k,T)/J2] 


P 

Bodner  and  Partom  further  expressed  as: 


D?  =  D2exp  -(Z2)n  n+1 

L  n  J 


(2.15) 


(2.16) 


This  expression  is  based  on  extensive  experimental  data  and 
has  been  modified  to  fit  results  found  by  several  researchers. 
Dq  is  the  limiting  value  of  the  plastic  strain  rate  in  shear. 
The  parameter  n  controls  strain  rate  sensitivity.  Z  is  the 
measure  of  material  hardening,  and  is  a  function  of  plastic 
work. 

Analysis  of  the  variables  in  the  equation  for  the  second 
invariant  of  the  plastic  strain  rate  show  that  Dq  is  an 
arbitrarily  chosen  value.  For  all  materials,  it  is  generally 

4 

chosen  as  Do=10  sec  ,  unless  very  high  rates  of  straining 
are  present.  The  values  of  n  chosen  affect  the  inelastic 
level  of  the  stress-strain  curves „  The  parameter  n  influ¬ 
ences  the  slope  of  the  stress-strain  rate  curve  and  is  there¬ 
fore  a  measure  of  strain  rate  sensitivity  [9] . 

Bodner’ s  elastic-viscoplastic  theory  is  based  on  an  in¬ 
ternal  material  state  variable  Z.  Bodner' s  Z  hardness 
parameter  is  a  basic  material  property  and  is  deformation 


history  dependent.  The  following  relationship  was  introduced 
for  the  Z  parameter: 


Z=Z(Wp)=Z1-(Z1-Z0)exp 


(2.17) 


The  variables  of  interest  in  the  above  relationship  include 

W  , Z.,Z  ,  and  m.  W  is  the  relative  amount  of  plastic  work 
p  x  u  p 

done  from  some  initial  state  and  is  mathematically  defined  at 
low  temperatures  to  be: 


V  \Sijfiijdt  (2-18)  • 

Z^  is  the  maximum  value  of  hardness,  ZQ  is  the  value  of  Z 

where  W  =0  (i.e  the  initial  state  point),  and  m  is  a  ma- 
P 

terial  constant,  hardening  rate  exponent  [10]. 

The  analysis  for  this  study  is  performed  on  a  material 
under  high  temperature  conditions.  It  is  therefore  necessary 
to  include  the  thermal  recovery  of  hardening  due  to  high  temp' 
erature.  The  plastic  work  oriented  term,  Wp  is  redefined 
as  follows: 


$ 


e_.?  dt 


Zrec  dt 

STz^zT 


(2.19) 


where 


(2.20) 


is  the  stable,  (i.e.,  non-work  hardened)  value  of 

at  a  given  temperature.  A  and  r  are  material  constants 

picked  to  match  creep  test  data.  The  thermal  recovery  term 

I [Z of  Eq  2.20  makes  a  negative  contribution  to  W  due 
rec  p 

to  the  negative  sign  on  A,  since  Z  is  always  greater  than 
or  equal  to  Z^. 

The  Bodner  constants  for  IN-100  at  1350  F  are  given  in 


Table  2.1. 


TABLE  2.1 

BODNER  COEFFICIENTS  FOR  IN-lOO  AT  1350°F 


Material 

Parameter 

Description 

Elastic  modulus 

Strain  rate  exponent 
Limiting  value  of  strain  rate 

Limiting  value  of  hardness 

Maximum  value  of  hardness 

Minimum  value  of  hardness 

Hardening  rate  exponent 

Hardening  recovery  coefficient 
Hardening  recovery  exponent 
(1  KBAR  *  100  MPa  «  14.504KSI) 


26.3x10  KSI 
(18.133xl0%>a) 


104  sec 


915.0KSI  (6304  MPa) 

1015. OKSI  (6993  MPa) 

600. OKS I  (4134  MPa) 

2.57KSr1  (.37273 
MPa-1) 

1.9xl0’3  sec'1 


Ill  Methods  of  Analysis 


The  Computer  Proaram 


A  two  dimensional  plane  strain  finite  element  program  named 


VISCO  was  used  throughout  this  study.  VISCO  accounts  for  both 
nonlinear  viscoplastic  material  behavior  and  changing  boundary 
conditions  due  to  crack  closure.  The  accuracy  of  the  program 
has  been  verified  by  Smail  [11]  and  Keck  [12] . 


The  Bodner-Partom  viscoplastic  constitutive  equations  in 
VISCO  are  solved  using  the  Gauss-Seidel  iterative  equation 
solver  with  overrelaxation,  eliminating  costly  stiffness 
matrix  factorization.  Time  integration  of  the  Bodner  equations 
is  accomplished  for  each  element  using  an  Euler  extrapolation 
scheme.  A  variable  time  step  algorithm  is  included  that  maxi¬ 
mizes  the  time  step  size  during  the  analysis  while  maintaining 
good  accuracy.  During  each  time  step,  equilibrium  tolerances 
are  checked.  If  the  tolerances  are  exceeded,  the  time  step  is 
reduced  until  equilibrium  is  obtained. 


To  completely  understand  the  material  behavior  of  IN-100, 
it  is  necessary  to  conduct  analysis  over  a  range  of  loads, 
frequencies  and  R-ratios.  The  data  can  then  be  applied  to  any 
life  duty  cycle  specified.  For  instance,  in  a  jet  engine 
turbine  disk,  there  can  be  frequent  load  cycles  of  varying 
frequencies  and  load  levels.  In  any  application,  there  is 
also  the  likelihood  of  components  being  subjected  to  compres¬ 
sive  loading. 

The  most  convienient  load  spectrum  for  engine  and  other 
applications  is  saw-toothed  load-time  pattern  with  constant 


VISCO' s  loading  function  was  modified  for  the  negative 
cycle  response  shown  in  the  lower  part  of  Fig  3.1  (see  Appendix 
A) .  Using  a  linear  load  equation  solver  to  model  the  load 
spectra  shown  in  Fig  3.1,  a  percentage  of  total  load  was  calc¬ 
ulated  at  each  time  step.  When  the  load  percentage  was  less 
than  zero,  the  load  is  applied  to  the  lower  edge  of  the  hole 
for  the  compact  tension  specimen  considered  in  the  study. 

When  the  load  percentage  is  positive,  the  load  is  applied  to 
the  top  edge  of  the  hole  in  the  compact  tension  specimen  (see 
Fig  3.4).  Node  displacement  and  element  strain,  stress  and 
the  Bodner  Z-hardness  were  output  at  requested  time  intervals. 

VISCO  was  also  modified  to  accomplish  crack  closure  and 
crack  opening.  The  boundary  condition  on  the  nodes  which  lie 
along  the  crack  edge  when  the  node's  displacement  is  positive 
is  such  that  there  is  zero  force  on  the  node.  When  the  node's 
displacement  is  calculated  to  be  negative  (indicating  closure) , 
the  boundary  condition  changes  to  one  of  zero  displacement, 
and  the  displacement  of  the  node  is  reset  to  zero. 


Finite  Element  Modelin 


Several  finite  element  models  were  utilized  throughout 
the  course  of  this  study.  The  first  was  the  2-element  uni¬ 
axial  model  shown  in  the  top  of  Fig  3.2.  This  model  was  used 
to  validate  the  changes  in  VISCO  and  to  provide  uniaxial  plane 
stress  solutions  for  the  frequencies  and  loads  shown  in  Table 


3.1.  Variations  in  the  computer  program  stress  and  strain 


increment  tolerances  were  tested  with  the  2-element  model,  and 
it  was  found  that  tighter  tolerances  than  those  used  in  the 
past  [12],  must  be  used  to  truly  characterize  the  material 
behavior  during  negative  load  cycling. 

TABLE  3.1 

UNIAXIAL  TEST  LOADS  AND  FREQUENCIES 

Test  Frequency  (Hz)  Max  Stresses  (ksi) 

2.5  130,  150,  165,  180,  200,  220 

0.167 

0.03 

The  next  model  used  was  the  10-element  uniaxial  model 
shown  in  the  lower  half  of  Fig  3.2.  This  model  was  used  to 
validate  the  accuracy  of  the  2-element  model  and  was  run  at 
a  frequency  of  0.167  Hz  and  a  maximum  stress  level  of  200  ksi. 

Two  dimensional  fatigue  modeling  was  accomplished  using  a 
standard  compact  tension  specimen  geometry  as  shown  in  Fig 
3.3.  Due  to  symmetry,  only  half  of  the  compact  tension  spec¬ 
imen  was  modeled  using  constant  strain  triangular  elements. 

The  finite  element  mesh  shown  in  Fig's  3.4  and  3.5  was  selected 
to  model  the  top  half  of  the  compact  tension  specimen.  This 
pattern  allows  for  unlimited  size  reduction  and  insures  that 
no  two  adjacent  elements  differ  in  size  by  more  than  a  factor 
of  2.  Except  for  elements  near  the  loading  pin  holes,  element 
aspect  ratios  varied  from  .5  to  1.0.  Elements  near  the  crack 
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TABLE  3.2 


COMPACT  TENSION  SPECIMEN  LOADS  AND  STRESS  INTENSITY  FACTORS  (K) 
(ksi  in)  Load  (lbs)  R-Ratio 


35 

1762 

0.1 

35 

1762 

-1.0 

45 

2265 

o 

0 

H 

1 

The  compact  tension  loads  are  shown  graphically  in  Fig  3.6. 

Verification  of  Computer  Program 

Validation  of  the  negative  cyclic  subroutine  and  research 
into  the  non-linear  response  on  IN-100  to  fatigue  loading  was 
accomplished  primarily  with  the  2-element  uniaxial  model. 

The  crack  closure  algorithm  was  verified  with  a  2-dimensional 
mesh  composed  of  543  elements  and  327  nodes.  This  mesh  was 
used  by  Keck  [12]  in  his  research,  which  considered  an  R-ratio 
of  0.1  and  is  of  known  accuracy.  Results  from  this  mesh, 
which  will  be  referred  to  as  mesh2A,  were  used  as  a  baseline 
in  the  effort  to  increase  the  number  of  load  cycles  permitted 
within  specified  computer  central  processor  time. 

Using  mesh2A  and  given  2,000  seconds  central  processor 
time  on  the  CDC  6600  computer,  only  4  complete  load  cycles 
were  possible.  A  new  mesh  was  created  which  is  composed  of 
382  elements  and  235  nodes.  This  mesh  is  herein  referred  to 
as  mesh382.  Using  mesh382  and  given  the  same  2,000  seconds 


central  processor  time,  13.5  load  cycles  are  possible.  The 
differences  between  mesh2A  and  mesh382  are,  the  number  of 
elements  (543  vs  382) ,  the  number  of  nodes  (327  vs  235) ,  and 
the  distance  ahead  of  the  crack  tip  where  the  element  area 
increases  (.078  in  vs  .048  in).  Comparison  of  the  results 
from  mesh382  with  the  results  from  mesh2A  show  less  than  3% 
error.  This  will  be  discussed  in  the  results  section  of  this 
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IV  Results  and  Discussion 

Results  and  comparisons  made  in  this  chapter  are  divided 
into  three  major  groups.  The  first  section  deals  with  extend¬ 
ing  existing  compact  tension  analysis  [12],  which  was  per¬ 
formed  with  an  R-ratio  of  0.1,  to  a  larger  number  of  load 
cycles.  The  results  were  highly  successful,  with  an  increase 
of  over  three  times  the  past  number  of  load  cycles  now  able 
to  be  analyzed.  The  second  section  is  concerned  with  the 
uniaxial  finite  element  specimen.  Here  is  the  first  attempt 
to  use  VISCO  for  negative  R-ratio  load  cycling.  The  program 
subroutines  are  verified  for  negative  loads,  and  the  accuracy 
of  the  solution  is  examined.  The  R-ratio  used  throughout 
this  portion  is  -1.0  and  various  cyclic  frequencies  and  load 
levels  are  examined.  The  final  section  deals  with  the  appli¬ 
cation  of  reversed  loading  to  the  compact  tension  specimen. 

The  R-ratio  used  in  this  section  is  -1.0,  the  load  frequency 
is  held  constant  at  2.5  Hz  and  two  load  levels  are  examined. 
The  crack  opening  displacements  and  residual  strain  illus¬ 
trations  yield  results  of  particular  interest. 

Positive  R-ratio  Compact  Tension  Specimen  Tests 

Results  and  illustrations  in  this  section  are  derived 
from  compact  tension  specimen  load  cycling  at  a  frequency  of 
2.5  Hz,  an  R-ratio  of  0.1,  and  a  maximum  stress  intensity  of 
K^=35  KSI  / in  .  The  discussion  is  broken  down  into  the 
following  subsections: 


vjv.vjv; 
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a)  Comparison  of  results  between  Mesh2A  and  Mesh382 

b)  The  amount  of  plastic  strain  accumulated  for  each  cycle 

c)  Vertical  (Y)  displacement  of  the  specimen  behind  the 
crack  tip 

d)  Cyclic  stress  strain  behavior  at  the  crack  tip  for  13 
load  cycles 

e)  Profiles  of  Z  hardness,  Y  stress  and  Y  strain  as  a 
function  of  distance  ahead  of  the  crack  tip  for  13  load  cycles 

Validation  of  the  accuracy  of  Mesh382  for  modeling  compact 
tension  specimen  behavior  was  accomplished  by  comparing  results 
obtained  with  results  given  from  Mesh2A,  which  is  known  to  be 
accurate.  Plastic  zone  sizes  after  2.5  cycles  are  compared  in 
Fig  4.1.  The  error  between  the  areas  of  the  two  plastic  zones 
is  less  than  2%,  which  is  considered  minor.  The  stress  fields 
ahead  of  the  crack  tip  for  the  two  meshes  after  2.5  cycles  are 
compared  in  Fig  4.2.  Again  the  differences  in  the  results 
between  Mesh382  and  Mesh2A  are  insignificant.  The  third  and 
final  check  c.i  the  accuracy  of  Mesh382  illustrated  in  this 
study,  is  a  comparison  of  the  specimen  corner  displacement 
through  the  second  load  cycle.  From  Fig  4.3,  it  is  obvious 
that  the  displacement  function  for  the  -  rner  of  the  compact 
tension  specimen  is  the  same  regardless  -f  which  finite  element 
mesh  is  used.  For  reasons  of  the  above  examples  and  other 
comparisons  of  strain,  Z  hardness  and  crack  opening  charac¬ 
teristics,  it  is  evident  that  Mesh382  will  accurately  model 
the  behavior  of  the  compact  tension  specimen. 
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BOTH  MESHES  SAME  POINT 


TANCE  AHEAD  OF  CRACK  TIP  (INCHES) 


o 


The  benefit  derived  from  using  Mesh382  for  modeling  the 
compact  tension  specimen  is  enormous.  Current  research  has 
been  able  to  achieve  only  a  moderate  number  of  load  cycles 
before  computer  operating  time  becomes  prohibitive.  For 
example,  Mesh2A,  which  has  543  elements  (327  nodes) ,  allows 
somewhat  less  than  four  complete  load  cycles.  With  equal 
accuracy,  Mesh382  is  able  to  complete  13.5  cycles  in  the  same 
amount  of  computer  time. 

The  first  area  to  be  examined  over  the  13  load  cycles, 
with  a  value  of  35  KSI  / in  ,  was  the  total  strain  over 
time  ahead  of  the  crack  tip.  Earlier  work  on  a  low  number  of 
load  cycles  indicated  that  a  stable  stress/strain  condition 
exists  at  the  crack  tip  after  the  first  load  cycle  (12] . 

There  appeared  to  be  a  constant  accumulation  of  plastic  strain 
at  each  cycle  in  the  material  ahead  of  the  crack  tip.  This 
behavior  has  little  physical  appeal,  because  materials  do  not 
normally  strain  without  bound  under  the  conditions  imposed  by 
the  compact  tension  specimen  test.  The  evidence  through  13 
load  cycles  refutes  the  concept  of  a  stable  stress/strain 
condition  ahead  of  the  crack  tip.  The  amount  of  plastic  strain 
accumulating  ahead  of  the  crack  tip  is  decreasing  in  an  expon¬ 
ential  fashion.  This  suggests  that  after  an  appropriate  number 
of  load  cycles,  the  total  strain  will  reach  a  stable  state  and 
no  further  plastic  straining  will  occur.  Table  4.1  illustrates 
the  value  of  the  total  strain  after  each  load  cycle  for  the 
first  element  ahead  of  the  crack  tip  along  with  the  change  in 
strain  from  the  previous  load  cycle.  This  data  is  presented 
in  Fig  4.4. 
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TABLE  4.1 


TOTAL  STRAIN  VALUES  AFTER  EACH  LOAD  CYCLE 

-3  -3 

Cycle  Strain (in/inxlO  )  Strain  Increase  (in/inxlO  ) 


0 

0.0 

0.0 

1 

14.03 

14.03 

2 

15.25 

1.22 

3 

16.07 

.827 

4 

16.58 

.51 

5 

17.0 

.42 

6 

17.33 

.33 

7 

17.60 

.27 

8 

17.84 

.24 

9 

18.09 

.25 

10 

18.29 

.20 

11 

18.49 

.20 

12 

18.67 

.18 

13 

18.81 

.14 

Extrapolation  of  the  data  in  Fig  4.4  indicates  that  after 
approximately  23  load  cycles,  no  more  plastic  straining  will 
take  place  (see  Fig  4.5).  Data  retrieved  at  this  point  could 
therefore  be  utilized  to  model  material  characteristics  during 
fatigue  tests  of  possibly  several  thousand  cycles. 

To  check  for  crack  closure,  vertical  (y)  displacements 
behind  the  crack  tip  were  monitored  during  the  unload  cycle 
for  each  of  the  13  load  cycles  examined  in  this  study.  As 
shown  in  Fig  4.6,  no  negative  displacements  are  observed. 


indicating  no  crack  closure.  This  is  to  be  expected  as  the 
minimum  load  still  applies  a  tension  in  the  specimen  near  the 
crack  tip  due  to  the  positive  R-ratio  of  0.1. 

The  displacement  of  the  corner  of  the  crack  mouth  for  the 
first  load  cycle  is  shown  in  Fig  4.7.  When  this  figure  is 
superimposed  over  the  crack  mouth  displacement  for  the  second 
load  cycle  shown  in  Fig  4.3,  it  is  found  that  the  slopes  of  the 
curves  are  identical.  Throughout  the  13  load  cycles  examined 
in  this  study,  there  was  no  change  in  the  plot  of  crack  mouth 
displacement  vs  load.  The  results  show  that  the  compact  ten¬ 
sion  specimen  behaves  elastically  at  points  on  the  specimen 
boundary  and  does  not  see  the  localized  inelastic  behavior  neer 
the  crack  tip. 

Fig  4.8  illustrates  the  stress/strain  behavior  for  the 
first  element  ahead  of  the  crack  tip.  Shown  are  the  first  four 
cycles,  and  then  cycle  numbers  8  and  13.  A  rapid  increase  in 
plastic  deformation  occurs  during  the  first  load  cycle.  Each 
successive  cycle  has  less  plastic  deformation  until,  as  shown 
earlier  in  Fig  4.5,  the  plastic  strain  per  cycle  goes  to  zero. 
After  about  23  cycles  then,  linear  elastic  stress  strain 
behavior  is  established.  The  constant  amplitude  of  stresses 
each  cycle  suggests  that  the  material  ahead  of  the  crack  tip 
is  essentially  in  a  stress  controlled  boundary  condition. 

The  behavior  of  the  Bodner  Z  hardness  material  parameter 
is  illustrated  in  Fig  4.9.  The  elements  immediately  ahead  of 
the  crack  tip  along  the  line  of  symmetry  reach  saturation 
value  within  the  first  load  cycle.  Elements  which  are  farther 
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out  than  2.4J5  of  the  crack  length  ahead  of  the  crack  tip  see 
no  change  from  the  initial  hardness  value  specified.  The 
elements  in  between  show  a  slight  increase  in  their  hardness 
value  through  approximately  three  load  cycles  where  no  more 
increase  is  indicated. 

The  stress  and  strain  fields  ahead  of  the  crack  tip  are 
shown  in  Figs  4.10  and  4.11.  The  stress  magnitude  is  rela¬ 
tively  constant  over  all  13  load  cycles,  supporting  the  sugges¬ 
tion  of  stress  controlled  boundary  condition.  The  strain  field 
shows  a  slight  increase  over  the  13  cycles  with  the  majority 
occurring  early  in  the  load  history.  This  is  again  showing 
the  trend  towards  a  stable  strain-time  history  after  approxi¬ 
mately  23  load  cycles. 
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PLASTIC  STRAIN  RATE  (IN/IN 


J'J  \9 


.008  .016  .024  .032  .040 

DISTANCE  AHEAD  OF  CRACK  TIP  I  INCHES  I 


DISTANCE  AHEAD  OF  CRACK  TIP  (INCHES) 


Negative  R-Ratio  Uniaxial  Results 

Three  major  areas  of  concern  were  examined  with  the  uni¬ 
axial  models  shown  in  Fig  3.2.  First,  the  subroutines  in  V.'.SCO 
for  reversed  loading  had  to  be  verified.  Next,  the  behavior 
of  IN-100  under  reversed  loads  was  examined.  This  data  was 
extracted  to  provide  a  baseline  for  the  behavior  of  the  2- 
dimensional  compact  tension  specimen.  Third,  and  last,  was  to 
examine  the  accuracy  of  VISCO  for  reversed  loading. 

An  accurate  stress/strain  solution  for  the  Bodner  Visco¬ 
plastic  equations  with  uniaxial  loading  under  high  stress 
levels  is  one  in  which  the  strain  increases  until  the  material 
fully  strain  hardens.  After  this  point,  the  stress/strain 
loops  for  each  cycle  will  overlay  each  other  exactly  [13,14]. 
The  accuracy  of  VISCO  is  partially  determined  by  the  input 
parameters  °tol  and  etol  which  are  the  stress  and  strain 
tolerances  respectively.  These  values  dictate  the  overall 
incremental  iteration  which  is  a  function  of  the  time  step. 

In  VISCO  the  time  step  dt,  for  time  step  i  is 


(4.1) 
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If  the  time  step  is  not  set  to  the  correct  value,  minor  error 
will  result  and  the  stress/strain  curves  will  not  exactly  over¬ 
lay  each  other.  The  stress  and  strain  tolerances,  unless 
otherwise  specified,  were  each  set  at  0.05. 

At  2.5  Hz  with  stress  levels  at  or  below  150  KSI  (see  Figs 
4.12,  4.13)  the  uniaxial  model  behaved  elastically.  A  steady 
state  stress/strain  solution  is  realized  immediately.  When 
the  stress  is  evaluated  to  165  KSI  (see  Fig  4.14),  some  small 
amount  of  plasticity  is  evident,  but  a  stable  stress/strain 
solution  is  present  immediately.  Fig  4.15  illustrates  the 
larger  plastic  strain  present  in  the  first  load  cycle  for  a 
stress  level  of  180  KSI.  The  material  hardens  so  that  after 
the  first  cycle,  a  stable  solution  is  seen  where  the  stress/ 
strain  loops  directly  overlay  each  other.  Some  variability  in 
the  solution  is  evident  for  stresses  of  200  and  220  KSI  (see 
Figs  4.16,  4.17).  This  is  due  to  the  slight  variations  in  the 
time  step  algorithm  discussed  in  equations  4.1  and  4.2. 
Neglecting  the  slight  variations,  a  stable  stress/strain  solu¬ 
tion  is  present  after  the  first  load  cycle. 

At  .167  Hz,  the  only  perfectly  elastic  stress/strain 
behavior  is  at  130  KSI  (see  Fig  4.18).  When  the  stress  level 
is  raised  to  150  KSI,  a  stable  solution  is  reached  after  the 
first  load  cycle,  as  illustrated  in  Fig  4.19.  Raising  the 
stress  level  to  165  or  180  KSI  introduces  some  minor  varia¬ 
bility  in  the  solution  as  shown  in  Figs  4.20  and  4.21.  It  is 
evident  that  the  material  should  see  a  steady  stress/strain 
solution  after  the  first  cycle.  Again  the  largest  plastic 


deformation  occurs  during  the  first  load  cycle.  At  stress 
levels  of  200  and  220  KSI  a  large  amount  of  plasticity,  and 
consequently  a  relatively  larger  amount  of  error  is  seen  in 
the  solution  (see  Figs  4.22  and  4.23).  According  to  the  Bodner 
Viscoplastic  equations,  the  stress/strain  loops  should  overlay 
one  another.  The  characteristics  of  the  largest  amount  of 
plasticity  occurring  during  the  first  cycle  is  evident  even 
in  these  two  figures. 

At  0.03  Hz,  stable  elastic-plastic  behavior  was  given  for 
stress  levels  up  to  and  including  165  KSI  (see  Figs  4.24-4.26). 
Little  variation  in  the  solution  is  noted  after  the  first 
cycle.  Input  stresses  of  180  and  200  KSI,  however,  produced 
a  non-stable  solution.  Figs  4.27  and  4.28  show  the  same  type 
of  increasing  strain  present  in  Figs  4.22  and  4.23.  Again, 
this  variation  is  due  to  the  time  step  algorithm  and  indirectly 
due  to  the  input  stress  and  strain  tolerances.  If  the  yield 
stress  of  IN-100  at  1350  F  is  set  at  130  KSI,  which  is  where 
the  stress/strain  curve  becomes  non-linear,  then  a  load  of 
220  KSI  is  significantly  higher  than  that  value,  and  if  applied 
for  more  than  a  small  amount  of  time,  should  lead  to  excessive 
strain  rates.  This  is  indeed  the  case  as  shown  in  Fig  4.29. 
Convergence  to  a  solution  was  not  possible  with  this  load 
input,  so  the  stress/strain  behavior  for  the  first  portion  of 
the  cycle  is  shown.  At  this  point,  the  time  step  became 
smaller  than  the  minimum  value  allowed  by  the  CDC  6600  computer 
and  no  further  computation  was  possible. 


Two  alternatives  were  examined  to  attempt  to  eliminate  the 
variability  in  the  solution  present  at  higher  stress  levels  and 
lower  frequencies.  First  the  10-element  model  of  Fig  3.2  was 
created  to  eliminate  any  errors  from  pin  and  roller  connections 
present  in  the  2-element  model,  which  are  not  realistic  for 
uniaxial  specimen.  The  stress  and  strains  of  elements  five 
and  six  were  averaged  and  then  plotted  over  four  load  cycles 
in  Fig  4.30.  The  input  stress  a^.^  frequency  were  identical  to 
that  of  Fig  4.22.  The  tightening  of  the  stress/strain  loops 
towards  the  point  where  they  would  begin  to  overlay  each  other 
indicates  that  the  10-element  model  is  a  step  towards  a  more 
accurate  solution.  The  other  alternative  was  to  decrease  the 
stress  and  strain  error  tolerances.  The  results  when  the 
tolerances  were  decreased  from  0.05  to  0.02  are  shown  in  Fig 
4.31.  Here  a  significant  move  towards  the  stress/strain  loops 
directly  overlaying  each  other  is  evident.  Since  the  most 
accurate  uniaxial  results  were  obtained  with  the  2.5  Hz  input 
frequency,  the  negative  cycling  of  the  2-dimensional  compact 
tension  specimen  was  accomplished  at  2.5  Hz. 

The  Bodner  Z  hardness  parameter  is  illustrated  in  Figs 
4.32-4.34.  The  value  of  Z  is  dependent  on  stress  level,  fre¬ 
quency  and  time.  In  Fig  4.32,  which  shows  the  Z  hardness 
values  for  2.5  Hz,  the  220  KSI  stress  level  saturates  within 
the  first  load  cycle  while  the  130  KSI  stress  level  remains  at 
the  input  value.  Similar  results  are  seen  in  Fig  4.33  for  a 
frequency  of  .167  Hz.  For  0.03  Hz,  shown  in  Fig  4.34,  the 
value  of  Z  reaches  its  maximum  within  the  first  load  cycle 
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Fig  4.16  Stress/Strain  Behavior  Over  Five  Cvcles 
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for  180  and  220  KSI  stress.  There  is  a  small  amount  of  de¬ 


crease  in  the  value  of  Z  during  the  elastic  portion  of  the  load 
cycle,  but  this  is  insignificant.  The  increase  in  Z  is  not 
instantaneous  under  any  of  the  load  conditions.  This  is  seen 
as  a  slope  in  the  step  function  in  Figs  4.32-4.34,  and  is 
characteristic  of  the  viscoplastic  action. 


Negative  R-Ratio  Compact  Tension  Specimen  Results 

Results  and  discussion  in  this  section  are  derived  from 
compact  tension  specimen  load  cycling  at  a  frequency  of  2.5  Hz, 
an  R-ratio  of  -1.0  and  maximum  stress  intensity  values  of  35 
and  45  KSlVin.  The  conclusions  reached  are  broken  down  into 
three  major  areas: 

a)  Crack  opening  displacements  behind  crack  tip 

b)  Profiles  of  the  stress  and  strain  fields  ahead  of  the 
crack  tip 

c)  Plastic  zone  size  estimations 

The  shape  of  the  crack  edge  as  the  load  goes  from  its  full 
tensile  value  to  its  full  compressive  value  during  the  second 
cycle  is  illustrated  in  Figs  4.35  and  4.36.  The  Y  displacement 
as  a  function  of  distance  behind  the  crack  is  identical  to  that 
of  a  positive  R-ratio  for  the  positive  portion  of  the  load 
cycle.  As  the  load  decreases  through  zero  and  becomes  compres¬ 
sive,  the  crack  begins  to  close.  The  elements  furthest  away 
from  the  crack  tip  close  first,  and  as  the  load  increases 


Fig  4.33  Z  Hardness  Over  Six  Cycles  at  0.167  Hz 
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in  compression,  the  crack  closes  closer  to  the  crack  tip. 

This  sequence  of  crack  closure  is  due  to  the  residual  plastic 
strain  around  the  crack  tip  providing  a  fulcrum  which  resists 
the  negative  displacement.  The  edge  of  the  crack  can  be 
thought  of  as  a  cantilever  beam  with  a  positive  angular  dis¬ 
placement  from  the  x-axis.  As  a  negative  force  is  applied  to 
the  end  of  the  beam,  it  will  go  to  zero  y-displacement .  Near 
the  fixed  support,  there  will  be  some  residual  positive 
y-displacement.  This  residual  y-displacement  provides  a 
radius  of  curvature  at  the  crack  tip  so  that  it  remains 
blunted,  even  under  totally  reversed  loads. 

The  stress  field  ahead  of  the  crack  tip  at  the  full  tensile 
load  value  is  shown  in  Figs  4.37  and  4.38.  This  profile  is 
steady  over  several  load  cycles.  This  type  of  stress  field 
has  been  discussed  by  several  of  the  authors  referenced  in 
this  thesis  [12,  13,  14].  Notice  that  for  K^=  45  KSI  / in 
(Fig  4.38),  the  stress  field  is  significantly  higher  than  for 
K1=  35  KSI  / in  .  A  larger  plastic  zone  is  established  as 
will  be  shown.  When  the  load  reaches  its  fully  compressed 
value,  the  stress  fields  shown  in  Figs  4.39  and  4.40  are 
established.  The  stresses  go  from  negative  to  positive  and 
back  to  negative  as  the  distance  ahead  of  the  crack  tip  in¬ 
creases.  This  profile  is  due  to  the  residual  positive  plastic 
strains  which  exist  even  at  the  fully  reversed  load  level. 

Note  that  the  magnitude  of  the  stresses  are  much  smaller  than 
the  corresponding  stresses  under  a  tensile  load.  Note  also 
that  the  reversed  load  for  K  =  35  KSI  /in  produces  a 


compressive  stress  less  than  that  for  a  K^=45KSI  /in  near  the 
crack  tip,  yet  just  the  opposite  is  occuring  away  from  the  tip. 

The  strain  field  ahead  of  the  crack  tip  under  a  tensile 
load  increases  in  magnitude  during  the  load  cycling  (see  Figs 
4.41  and  4.42).  It  is  obvious  that  the  larger  value  pro¬ 

duces  a  strain  function  of  greater  magnitude.  A  fairly  con¬ 
stant  rate  of  strain  increase  per  cycle  is  evident  for  both 

values.  This  indicates  that  the  stable  strain  solution 
for  reversed  loading  may  be  at  a  larger  number  of  cycles  than 
for  positive  R-ratio  load  spectra.  When  the  load  is  at  its 
full  negative  value  the  strain  immediately  ahead  of  the  crack 
tip  is  initially  negative  for  K^=  35  KSI  /in  ,  but  becomes 
positive  within  four  cycles  (see  Fig  4.43).  The  strain  com¬ 
ponent  for  K^=  45  KSI  /in”  is  always  positive  immediately 
ahead  of  the  crack  tip.  These  residual  positive  strains  are 
due  to  the  large  plastic  overload  which  occurs  during  the  first 
load  cycle.  The  stress  profile  ahead  of  the  crack  tip  under 
compressive  loading  is  governed  by  these  residual  strains 
interfacing  with  the  negative  input  loads. 

The  effect  of  load  magnitude  on  the  size  of  the  plastic 
zone  is  examined  in  Figs  4.45  and  4.46.  These  two  figures 
illustrate  the  approximate  shape  and  size  of  the  plastic  zone 
when  the  load  is  at  its  maximum  tensile  value.  An  element  is 
assumed  to  have  undergone  plastic  strain  if  the  magnitude 

is  greater  than  1x10  ^  in/in.  The  size  of  the  plastic  zone  is 
very  stable  over  load  cycling.  The  slight  increases  in  strain 
shown  in  Fig  4.41  and  4.42  are  not  large  enough  to  drive  any 
of  the  elements  adjacent  to  the  plastic  zone  into  plasticity 
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over  the  number  of  load  cycles  examined.  When  the  load  is 
reversed,  the  plastic  zone  is  compressed  into  an  elliptically 
shaped  region  (see  Figs  4.47  and  4.48).  The  area  of  this 
region  is  approximately  one-twelfth  of  the  area  of  the  plastic 
zone  under  a  tensile  load.  This  region  of  significant  posi¬ 
tive  y  strain  which  exists  under  the  compressive  load  is  due 
to  the  large  plastic  overload  which  occurs  during  the  first 
load  cycle.  The  type  of  material  behavior  which  occurs  under 
this  condition  could  be  expected  to  follow  the  same  pattern 
as  that  which  is  seen  when  a  specimen  under  positive  loading 
undergoes  an  application  of  an  overload.  The  material  sur¬ 
rounding  the  plastic  zone  applies  a  compressive  force  which 
decreases  its  size.  The  resulting  residual  compressive  stresses 
ahead  of  the  crack  tip  will  aid  in  retarding  crack  growth. 
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V  Conclusions 


The  results  of  the  computations  carried  out  using  the 
Bodner  viscoplastic  constitutive  equations  to  model  the  be¬ 
havior  of  IN-100  at  elevated  temperatures  yield  the  following 
conclusions  for  positive  R- ratios. 

1 .  The  increase  in  strain  realized  after  each  load  cycle  under 
an  R-ratio  of  .1  decreases  in  a  manner  which  indicates  that 
after  23  load  cycles,  the  material  will  no  longer  undergo  any 
more  plastic  straining.  All  data  retrieved  at  this  point  would 
therefore,  be  applicable  to  fatigue  tests  of  several  thousand 
cycles. 

2.  The  large  majority  of  plastic  straining  occurs  within  the 
first  three  load  cycles  regardless  of  the  R-ratio  utilized. 

3.  The  stress  field  ahead  of  the  crack  tip  remains  relatively 
constant  after  one  to  three  load  cycles  regardless  of  the 
R-ratio  utilized. 

4.  The  size  of  the  plastic  zone  is  load  dependent. 

5.  The  cyclic  behavior  of  the  compact  tension  specimen  near 
the  crack  tip  is  stress  controlled. 

6.  The  far  field  overall  behavior  of  the  compact  tension 
specimen  is  elastic  and  independent  of  load  level. 

The  following  conclusions  can  be  made  based  on  the  negative 
cyclic  loading  of  a  compact  tension  specimen  after  observing 
the  results  for  five  load  cycles: 

1.  The  edges  of  the  crack  do  not  meet  when  the  load  is  released 
due  to  the  large  plastic  overstress  in  the  first  load  cycla. 


2.  The  edges  of  the  crack  come  into  contact  when  a  compressive 
load  is  applied,  but  an  open  region  remains  immediately  behind 
the  crack  tip  due  to  the  large  plastic  overstress  in  the  first 
cycle. 

3.  An  elliptically  shaped  region  of  positive  Y-strain  is 
present  in  the  vicinity  of  the  crack  tip  even  under  totally 
reversed  loads.  The  compressive  load  decreases  the  area  of 
the  significant  positive  Y-strain  by  a  factor  of  twelve. 

4.  The  crack  tip  remains  blunted  under  reversed  loading  which 
other  researches  indicate  leads  to  a  higher  fatigue  life. 
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APPENDIX  A 


The  computer  program  VISCO  was  modified  for  negative  cyclic 
loading  in  several  areas.  The  first  change  was  made  in  the 
"read  input  data"  section.  Previously  loads  were  only  applied 
to  the  top  of  the  hole  when  the  input  load  is  applied.  Since 
the  load  was  now  to  be  applied  at  either  the  top  or  the  bottom 
of  the  hole  (depending  on  the  algebraic  sign  of  the  load) ,  the 
node  numbers  along  the  bottom  of  the  hole  needed  to  be  included 
as  part  of  the  input  data.  These  nodes  are  identified  by  the 
variable  name  NFN(I)f  where  I  goes  from  one  to  the  number  of 
nodes  to  be  loaded  in  compression.  The  numbers  of  the  nega¬ 
tively  loaded  nodes  are  read  in  line  980. 

The  load  is  shifted  from  the  top  of  the  hole  to  the  bottom 
of  the  hole  and  back  to  the  top  in  the  "LOAD"  subroutine. 

This  shift  of  load  depends  on  the  sign  of  the  input  load.  If 
the  load  is  positive,  the  force  is  applied  to  nodes  NFA(I) , 
which  are  along  the  top  of  the  hole.  If  the  load  is  negative, 
the  force  is  applied  to  nodes  NFN(I),  which  are  along  the 
bottom  of  the  hole.  The  algorithm  for  shifting  the  load  is 
contained  in  lines  6470-6550. 

To  accomplish  crack  closure  and  crack  opening,  changes 
were  made  in  the  "SOLVE”  subroutine.  Nodes  along  the  crack 
edge  (91-97)  could  not  be  allowed  to  have  negative  displace¬ 
ments  and  must  also  have  a  zero  force  boundary  condition  for 
positive  loads.  A  test  was  incorporated  for  the  above  nodes 
to  determine  if  the  displacement  was  negative,  and  if  so,  set 


it  equal  to  zero.  The  algorithm  for  crack  closure  and  opening 
is  contained  in  lines  7060-7090. 

The  remaining  alterations  were  done  in  the  "CYCLIC"  sub¬ 
routine.  The  changes  allow  any  R-ratio  and  frequency  to  be 
modeled.  TTOP  and  TBOT  are  variables  to  indicate  at  what  time 
the  load  should  be  at  its  maximum  and  minimum  values  respec¬ 
tively.  The  load  percentage  PP  is  then  calculated  considering 
the  point  of  time  within  the  cycle  given  by  TLOC  and  a  linear 
load  and  unload  profile.  Changes  in  subroutine  CYCLIC  are  in 
lines  12130-12200. 
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studies  have  incorporated  compressive  loads  or  analyzed  the 
fatigue  effect  of  a  negative  R-ratio  on  the  fatigue  characteris¬ 
tics  of  the  superalloy  in  IN-100. 
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This  study  involves  extending  existing  analysis  of  the  stress 
field  and  plastic  zone  ahead  of  a  crack  tip  in  a  compact  tension 
specimen/  through  a  larger  number  of  load  cycles  and  examining 
the  nearly  unexplored  area  of  compressive  loading  of  a  crack  in 
a  superalloy.  A  USAF  Materials  Laboratory  finite  element  com¬ 
puter  program  named  VISCO  was  used  for  this  study.  The  Bodner- 
Partom  viscoplastic  constitutive  equations  for  describing  the 
material  behavior  were  utilized.  Load  spectra  included  various 
frequencies  with  R-ratios  of  0.1  and  -1.0  (zero  mean  load). 
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